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ABSTRACT 

We quantify the fraction of galaxies at moderate redshifts (0.1 < z < 0.5) that appear red-and-dead 
in the optical, but in fact contain obscured star formation detectable in the infrared (IR), with the 
PRIsm MUlti-object Survey (PRIMUS). PRIMUS has measured ~ 120,000 robust redshifts with a 
precision of er 2 /(l + z) ~ 0.5% over 9.1 deg 2 of the sky to the depth of i ~ 23 (AB), up to redshift 
z ~ 1. We specifically targeted 6.7 deg 2 fields with existing deep IR imaging from the Spitzer Space 
Telescope from the SWIRE and S-COSMOS surveys. We select in these fields an i band flux-limited 
sample (i < 20 mag in the SWIRE fields and i < 21 mag in the S-COSMOS field) of 3310 red-sequence 
galaxies at 0.1 < z < 0.5 for which we can reliably classify obscured star-forming and quiescent galaxies 
using IR color. Our sample constitutes the largest galaxy sample at intermediate redshift to study 
obscured star formation on the red sequence, and we present the first quantitative analysis of the 
fraction of obscured star-forming galaxies as a function of luminosity. We find that on average, at 
L ~ L* , about 15% of red-sequence galaxies have IR colors consistent with star- forming galaxies. The 
percentage of obscured star-forming galaxies increases by ~ 8% per mag with decreasing luminosity 
from the highest luminosities to L ~ 0.2L*. Our results suggest that a significant fraction of red- 
sequence galaxies have ongoing star formation and that galaxy evolution studies based on optical color 
therefore need to account for this complication. 

Subject headings: galaxies: fundamental parameters (classification, colors, luminosities, masses, radii, 
etc.) — galaxies: distances and redshifts - galaxies: evolution - stars: formation - 
dust, extinction 



1. INTRODUCTION 

Our understanding of the formation and evolution 
of galaxies has improved tremendously the last few 
decades as a result of multiple redshift surveys (e.g ., 
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surveys have demonstrated 
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that the distribution of galaxy populations is bimodal 
the opt ical color-magnitude space up to 
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iStrateya et al.1 120011: iBlanton et al.l 12003a; 
Baldrv et al.l 120041; iBalogh et alj|2004t jBell et at|[2 004b; 
Wilhner et all 120061; jWuvts et al.l 120071: iWilliams et al I 
2009i; IBrammer et al.l l2009. among others) . The underly- 
ing physical reason is the short lifetime of massive stars. 
Once star formation stops in a galaxy, hot massive stars 
quickly die, the integrated spectral energy distribution 
(SED) of the galaxy becomes dominated by cool giants, 
and its optical color turns from blue to red. The bimodal- 
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ity of blue cloud and red sequence in the color- magnitude 
diagram therefore approximately translates into that of 
star-forming (SF) and quiescent galaxies. 

Recent studies of the mass and luminosity functions 
of red-sequence galaxies at different redshifts show that 
the stellar mass density on the red sequence has in- 
creased by a factor of ~ 2 overa ll since z ~ 1 (e.g., 
iBell et allboTjIHlFaber et al.|[2007l ). This build-up could 
result from some combination of migration from the blue 
cloud (due to mergers, feedback or other mechanisms for 
turning off star formation) and gas-poor, "dr y" merging 
from smaller-m ass red-sequence galaxies fe .g.. IBell et al. | 
" 2004H, l2006allbT). Other studies (e.g., iCimatti et all 
20061: Brown et all [20071: IC00I et alj|200g ) find that the 



growth of the red-sequence populations depends strongly 
on luminosity: whereas many L* red-sequence galaxies 
assembled in the time since z ~ 1, at the massive end 
(L > 3L*) the mass and luminosity functions are consis- 
tent with passive evolution. They therefore argue that 
merging (dry or wet) does not dominate the evolution of 
massive red galaxies since z ~ 1. 

One complication that current studies of red-sequence 
galaxies ignore is that optically red galaxies can have in 
situ, but dust-obscured star formation. Because of this 
dust-obscuration, optical color can be a misleading proxy 
for star formation. Dust extinction reshapes the SED of 
a galaxy, absorbing photons at high energy (e.g., ultra- 
violet, UV) and re-emitting the energy at longer wave- 
lengths (e.g., infrared, IR). Therefore the red sequence 
consists of both truly quiescent but also obscured SF 
galaxies; it is extremely important to quantify the com- 
position of the red sequence to obtain an accurate un- 
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derstanding of galaxy evolution. 

There are several ways to investigate the contribution 
of obscured SF galaxies on the red sequence. Assuming 
inclination is we ll-correlated with d ust extin ction along 
the line-of-sight. IMaller et al.l (120091 see alsolShao et al.l 
120071: lUnterborn fc RvdeiJl2008t IBailin fe Harrisl l2008b) 
find the ratio of red-to-blue galaxies changes from 1:1 to 
1:2 when moving from observed to inclination-corrected 
color, for l ocal galaxies in the Sloan Digital Sky Sur- 
vey (SDSS, lYork et all [2000). Using deep imaging from 
the Hubble Space Telescope (HST) from the Galaxy 
Evolution from Morphology and SED project (GEMS, 
iRix et all I2004D a nd photometric r edshifts from the 
COMBO-17 survey (|Wolf et al.ll200l . lBell et al.l (|2004aD 
find that less than 13% of red-sequence galaxies are 
edge-on dusty disk galaxies at z ~ 0.7. One can also 
correct dust extinction by fitting multi-wavelength data 
or SEDs with stellar population synthesis models and 
allowing dust extinction to be a fre e paramete r . Us - 
ing this method on COMBO-17 data, iWolf et al.l (|200l 
find a rich component of dusty SF galaxies contaminat- 
ing the red sequence in the cluster Abell 901/902 at 
z ~ 0.17. The mid-IR (MIR) emission from dust of- 
fers another method for selecting obscured SF galaxies, 
though the emission can als o originate f r om st rong ac- 
tive galactic nuc l ei (AGNs). IBrand et all (120091 see also 
Coia et al.l 120051: iDavoodi et al.l 120061: iRodighiero et al.l 
20071: ISaintonge et al.l 120081) discover ~ 10% of red- 
sequence galaxies at 0.15 < z < 0.3 exhibit strong MIR 
24 /an emission with / 2 4 > 300 ^Jy. We can also use 
IR color to separate SF and quiescent galaxies. At wave- 
lengths < 10 /im, warm dust emission and strong poly- 
cyclic aromatic hydrocarbon (PAH) features dominate 
the SED of SF galaxies, producing much r edder IR col- 
ors th an expected from quiescent galaxies. IBrand et al.l 
(2009) use this method and find that, among the ~ 10% 
of red-sequence galaxies with / 2 4 > 300 /^Jy in their sam- 
ple, about 64% have IR colors consistent with SF galax- 
ies, with the rest consistent with quiescent galaxies (31%) 
and AGNs (5%). 

Quantitative studies of obscured star formation at in- 
termediate and high re dshifts are usually limited to small 
samples. For example, Wolf et ail (|2005l ) select 462 red- 
sequence galaxies in a specific environment, the Abell 
901/902 super-cluster system. Even though deep multi- 
wavelength data including space-based IR imaging are 
now available across a large area of the sky, the lack of 
corresponding large redshift surveys hampers the system- 
atic quantitative studies of the contribution of obscured 
SF galaxies on the red sequence. We have recently com- 
pleted such a redshift survey, the PRIsm MUlti-object 
Survey (PRIMUS0), to provide the redshifts required to 
fully realize the science potential of the existing deep 
multi-wavelength imaging data. 

PRIMUS has measured ~ 120,000 robust redshifts 
over 9.1 deg 2 of the sky to the depth of i ~ 23 (AB) up to 
redshift z ~ 1. We specifically targeted areas where there 
exist deep multi-wavelength imaging data, including IR 
data from the Spitzer Space Telescope (jWerner et al.l 
120041) . making our dataset ideal for the study of the ob- 
scured star formation on the red sequence. Using the IR 
data, in conjunction with the redshifts from PRIMUS, 



we here present the first quantitative study of the frac- 
tion of obscured SF galaxies on the red sequence as a 
function of luminosity at redshifts 0.1 < z < 0.5. 

This paper proceeds as follows. In Section^ we briefly 
introduce PRIMUS and describe the Spitzer data we 
use. In Section [31 we describe in detail the classification 
scheme for AGNs, star-forming galaxies, and quiescent 
galaxies. We present the results in Section [4] In Section 
El we discuss the composition of the red sequence. We 
summarize our principal conclusions in Section |6l 

Throughout this work, we adopt a ACDM cosmo- 
logical model with f2 m = 0.3, fl\ = 0.7 and Hq = 
lOO/i km s _1 Mpc -1 . If not specifically mentioned, all 
magnitudes are on the AB system. 

2. DATA 

2.1. PRIMUS 

PRIMUS is a unique low-resolution spectroscopic in- 
termediate redshift survey. We conducted the survey us- 
ing the IMACS instrument on the Magellan I (Baadc) 
Telescope at Las Campanas Observatory. With a prism 
and slitmasks, we obtained low-resolution spectra for 
~ 300, 000 objects to the depth of i ~ 23 out to red- 
shift z ~ 1 in 39 nights, covering 9.1 deg 2 of the sky. 
The low-resolution of the prism (R ~ 30 — 150) allows 
us to observe ~ 2500 objects in one single mask with 
a field of view 0.18 deg 2 . By carefully fitting the low- 
resolution spectra with galaxy spectral templates, we 
have obtained ~ 120,000 robust redshifts, to the pre- 
cision of cr z /(l + z) ~ 0.5%. PRIMUS therefore is the 
largest faint galaxy redshift survey performed to date. 

When designing the survey, we specifically targeted 
fields with existing deep multi-wavelength imaging, in- 
cluding optical imaging from various ground-based deep 
surveys, near-UV (NUV) and far-UV (FUV) from 
the G alaxy Evolution Explorer (GALEX, iMartin et al.l 
l2005f ). IR from the Spitzer Space Telescope, and X-ray 
from Chandra and/or XMM telescopes. 

We created homogeneous band-merged catalogs across 
all our fields. With the redshifts and the band-merged 
catalogs, we derived ab solute magnitudes usin g the 
kcorrect package (v4.1.4 lBlanton fc Roweisl[2007T) . Be- 
cause we focus on redshifts 0.1 < z < 0.5, here we use 
the absolute magnitude s at uqriz bands shi fted blue- 
ward to z = 0.3 (e.g., IBlanton et aT]|2003aI ). We de- 
note these bands with a preceding superscript 0.3, e.g., 
3 z. For a galaxy exactly at redshift 0.3 with apparent 
magnitude i, the ^-correction from i to Mo .3j is simply 
— 2.51og 10 (l + z) independent of galaxy SED. This filter 
choice therefore minimizes the error in the ^-corrections. 

When fitting for the redshift, we also fit the spectra 
with stellar and broad-line AGN spectral templates and 
classify objects into stars, AGNs, and galaxies based 
upon the x 2 °f the best fits. In ~ 5.5 deg of the 
9.1 deg 2 PRIMUS fields, we have X-ray data from Chan- 
dra and/or XMM telescopes. Here we select our sample 
only from objects that are classified as galaxies and do 
not have matched X-ray detections. 

For more details on the P RIMUS survey, we refer the 
reader to the survey papers (jCoil et al.ll2010L Cool et al., 
in prep.). 



http: / / cass.ucsd.edu/~acoil / primus / 



2.2. Spitzer Space Telescope 
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In - 7.2 deg 2 of the 9.1 deg 2 fields covered by PRIMUS 
we have IR im aging data from The Infrared Array 
Camera (TRAC. iFazio et al.l I2004D and The M ulti-band 
Imaging Photometer (MIPS. lRieke et al.ll2004f ) onboard 
Spitzer. IRAC is a four-channel instrument that pro- 
vides simultaneous broadband images at 3.6, 4.5, 5.8 and 
8.0 /im with unprecedented sensitivity. MIPS provides 
images at longer wavelengths, 24, 70 and 160 //m with 
5cr rms depths of 0.3 mJy, 25 mJy and 150 mJy, respec- 
tively. Here we use the ~ 6.7 deg 2 that has coverage from 
all four channels of IRAC and 24 /im band of MIPS. For 
- 0.5 deg 2 of the DEEP2 (|Davis et al.ll200l field with 
the IRAC imaging, only the 3.6 /im and 4.5 /im imag- 
ing is available; we omit this field, because we require 
all four channels for this paper. The data we use con- 
sist of two Legacy Programs of Spitzer, 5.7 deg 2 in three 
fields of Spitzer Wide-area Infrared Extragalactic Sur- 
vey (SWIREFI: ILonsdale et aT1l2003l) . and 1 deg 2 in the 
Spitzer Deep Survey o f the HST COSMOS 2-Degree ACS 
Field (S-COSMOljB; ISanders eFall 120071 ) . We use the 
data release 2 and 3 of the SWIRE survey, and the pub- 
lic catalog in the NASA/IPAC Infrared Science Archive 
for the S-COSMOS survey. In most of the SWIRE fields, 
to at least ha depths, the flux limits in the four IRAC 
channels are 10, 10, 43, and 40 /xJy (21.4, 21.4, 19.8 and 
19.9 AB ma£3). The flux limits in the S-COSMOS field 
are 0.9, 1.7, 11.3, and 14.6 /iJy (24.0, 23.3, 21.3 and 21.0 
AB mag), in the four IRAC channels, respectively. 

At wavelengths < 10 /im, warm dust emission and 
strong PAH features, e.g, 3.3, 6.2, 7.7 and 8.6 /xm, dom- 
inate the SEP of a typical SF galaxy (|Li fc Draindl2001l : 
iSmith et al.ll2007l ). Moving to higher redshift, the 6.2, 7.7 
and 8.6 /im PAH emission bands shift through the IRAC 
8.0 /im channel, producing much redder IRAC colors 
than expected from a typical quiescent early-type galaxy. 
Red-sequence galaxies whose star formation is obscured 
by dust have IRAC colors consistent with normal SF 
galaxies. Strong AGNs, on the other hand, occupy a 
different locus from regular SF and quiescent galaxies 
in the IRAC colo r -color diagram (lEisenhardt et al.ll2004j : 
ILacv et al.l l20f)i IStern et all [2005J). The IRAC colors 
therefore are extremely convenient and powerful for clas- 
sification of AGNs, SF galaxies, and quiescent galaxies. 
We discuss this in detail in Section [3J 

3. SOURCE SELECTION 
3.1. Flux-limited sample in i band 

In order to use IR color to identify SF galaxies, we re- 
strict our sample to the PRIMUS fields with IR coverage 
(6.7 deg 2 ). We also only select our sample from galaxies 
with the highest confidence (Q = 4, see Cool et al. in 
prep, for details on the confidence level Q). 

We also need to select appropriate flux limits in the i 
band and the 3.6 /im band, which are necessarily related 
to the depth of the coverage in the 8.0 /im band. The 

10 http:/ /swire. ipac.caltech.edu/swire/swire. html 

11 http://irsa.ipac.caltech.edu/data/SPITZER/S-COSMOS/; 
COSMOS represents the Cosmo-logical Evolution Survey and ACS 
is the Advanced Camera for Surveys onboard HST. 

12 We use AB magnitudes throughout this work but also 
give Vega magnitudes when necessary. The conversion factors 
(AB-Vega) for IRAC 3.6, 4.5, 5.8, and 8.0 pirn are 2.79, 3.26, 
3.73, and 4.40 mag, respectively. 
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Fig. 1. — Selection of the i band flux-limited sample. We present 
[3.6] vs. i for red-sequence galaxies. The detection limits in the 
8.0 /im channel are 19.9 (21.0) mag in the SWIRE (S-COSMOS) 
fields, effectively much brighter than the limits in the 3.6 /im chan- 
nel. Many 3.6 fim detections are not detected in the 8.0 /im chan- 
nel. We can only use the IR color [3.6] — [8.0] to perform re- 
liable star-forming and quiescent galaxy classification for objects 
with bright 3.6 /im fluxes with [3.6] < 20 (21) mag. We show 
here that most of the red-sequence galaxies with i < 20 (21) mag 
have [3.6] < 20 (21) mag. We therefore only select galaxies with i 
brighter than 20 (21) mag in the SWIRE (S-COSMOS) fields. In 
Section 13.41 we show that these cuts enable us to classify almost 
all non-detections in 8.0 /im channel as quiescent galaxies in this 
flux-limited red-sequence sample. 

flux limits used in the 8.0 /im channel are 19.9 mag in 
the SWIRE fields and 21.0 mag in the S-COSMOS field, 
effectively much brighter than those in the 3.6 /im chan- 
nel (21.4 and 24.0 mag). Thus, in order to place inter- 
esting limits on the specific star-formation rate (SSFR) 
using the [3.6] — [8.0] color, we must restrict our sam- 
ple to bright enough galaxies. In Section 13.41 we will 
show that galaxies with SSFR of ~ 10~ 10 yr _1 have 
colors [3.6] — [8.0] > and quiescent galaxies have 
[3.6] - [8.0] < 0. Our definition of "star-forming" will 
be loosely related to this SSFR, and this color cut de- 
fines how bright our sample needs to be. We thus require 
the upper limits of the [3.6] — [8.0] colors of the 8.0 /im 
non-detections to be < so that we can reliably classify 
them into quiescent galaxies, or in other words that the 
[3.6] magnitudes should be brighter than ~ 20 (21) mag 
in the SWIRE (S-COSMOS) fields. 

We achieve this by applying appropriate limits in the 
i band. In Figure [TJ we present the [3.6] versus i for the 
red-sequence galaxies (defined in the next section) , where 
we have shown [3.6] = 20 and 21 mag with two horizon- 
tal lines. We see that for most of the red-sequence galax- 
ies with i < 20 (21) mag, they have 3.6 /im magnitudes 
[3.6] < 20 (21). We therefore define a flux- limited sample 
with i < 20 and i < 21, in the SWIRE and S-COSMOS 
fields, respectively. We will further justify these cuts in 
Section 13.41 where we explicitly demonstrate that after 
these cuts, virtually all 8.0 /im non-detections are quies- 
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Fig. 2. — Color-magnitude diagram of galaxies with i < 22.5. We 
use the absolute magnitudes in g and i bands shifted blueward to 
2 = 0.3 ( 3 g and 3 j) to minimize the errors in X-corrections. We 
show the diagram in two redshift bins: 0.1 < z < 0.3 (lower panel) 
and 0.3 < z < 0.5 (upper panel). The red solid lines and the blue 
solid lines indicate the positions of the red sequence and the blue 
cloud. The two green dashed lines represent the selection criteria 
(Equation (1) and (2)) we adopt to select rcd-scqucncc and blue- 
cloud galaxies, at median redshift in each bin. Galaxies between 
them are identified as green-valley galaxies. At different redshifts, 
we fix the slope of the criteria but allow the intercept to evolve 
(redden) by 0.4 per unit redshift. 

cent galaxies. 

Beyond redshift 0.5, the 6.2 fim PAH feature shifts out 
of the 8.0 /j,m channel and we can not use the [3.6] — [8.0] 
color to reliably separate SF and quiescent galaxies. We 
also do not have many galaxies below redshift 0.1 in the 
PRIMUS fields due to the small volume. Therefore, we 
focus our analysis on the redshift range 0.1 < z < 0.5. 
We will come back to this in Section 13.41 In total, the 
flux-limited sample within 0.1 < z < 0.5 includes 7258 
galaxies. 

3.2. Red sequence, blue cloud, and green valley 

To define the red sequence, blue cloud, and green val- 
ley, we use all galaxies with i < 22.5 with secure red- 
shifts (Q = 4) in PRIMUS. PRIMUS targeted in i band 
unless there is no i band imaging available, in which case 
we used R band. We thus choose to use the magnitude 
Afo.3j — 51og 10 h, and the color Mo.3 g — Mo.3 i when defin- 
ing samples based on optical color. At each luminosity, 
we fit the color distribution with a two-gaussian model. 
We then linearly fit the two sets of peaks to obtain the 
slopes and intercepts of the red sequence and the blue 
cloud. To define the edges of green valley, we move the 
fitted line for the red sequence (blue cloud) blueward 
(redward) until it includes 60% more galaxies above (be- 
low) the line. This method allows us to conservatively 
select red-sequence and blue-cloud galaxies, without con- 
tamination from green-valley galaxies. Due to passive 
evolution, the colors of galaxies at lower redshift are red- 
der than those at higher redshift. Therefore, we keep 



Fig. 3. — AGN classification for detections in both the 5.8 fim and 
the 8.0 /an channels in the flux-limited red-sequence sample. We 
show the IRAC color-color diagram with [3.6] - [4.5] vs. [5.8] - [8.0] . 
T he dotted wedge sh ows the criteria for AGN selection suggested 
by Stern ct al. (2005). The dot-dashed line represents the colors of 
power-law spectra f(v) = v~~ a of AGNs with the index a ranging 
from 0.5 to 3.0. Galaxies that fall into the wedge are identified 
as AGNs. For non-AGN galaxies, we use the [3.6] — [8.0] color 
to classi fy th em into star-forming galaxies and quiescent galaxies 
(Section UJJl. We show the SF galaxies in blue filled triangles and 
the quiescent galaxies in red filled stars. 

the slope fixed throughout the redshift range, but allow 
the intercept to evolve linearly in redshift. By fitting 
the color-magnitude diagram at different redshifts, we 
measure the color reddening is ~ 0.4 per unit redshift in 

Mo.3g — Mo.3^. 

Finally, the criteria we adopt are as follows. We require 
a red-sequence galaxy to have the color: 

Mo.3 g — Mo.3j > 

1.73 - 0.11(Mo.3i - 51og 10 h + 20) 

-0.4(2-0.3), (1) 

and a blue-cloud galaxy to have the color: 

Mo.3g — Afo.3j < 

1.40 - 0.10(Mo.3 j - 5 log 10 h + 20) 

-0.4(z-0.3). (2) 

Between these two criteria, we identify objects as green- 
valley galaxies. As an example, in Figure [H we show the 
color-magnitude diagram of all galaxies with i < 22.5 at 
0.3 < z < 0.5 (upper panel) and at 0.1 < z < 0.3 (lower 
panel). The two solid lines in each panel represent the 
positions of the red sequence and the blue cloud, and the 
dashed lines indicate the edges of the green valley, i.e., 
the criteria defined by Equation (1) and (2) at median 
redshift in each bin. 

Finally, in the flux-limited sample, we have 3310 red- 
sequence galaxies, 2497 blue-cloud galaxies and 1451 
green- valley galaxies. 
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Fig. 4. — AGN classification for non-detections in the 
5.8 fim channel that are detected in the 8.0 (im channel in the flux- 
limited red-sequence sample. We use the flux limits of [5.8] < 19.8 
(21.0) mag in the SWIRE (S-COSMOS) fields to calculate the lower 
limits of their [5.8] — [8.0] colors. We can safely treat all the non- 
detections as non-AGNs. 

3.3. AGN identification 

For the purposes of this paper, we are primarily inter- 
ested in star-forming galaxies, and attempt to identify 
and exclude the few AGNs in the sample. To do so, 
we use a combination of the PRIMUS spectra, matching 
X-ray data, and IR colors of the galaxies. 

As stated in Section 12.11 we only select objects clas- 
sified as galaxies and did not include those classified as 
broad-line AGNs in the PRIMUS spectroscopy. For Ar- 
ray AGNs, we exclude galaxies with matched A-ray de- 
tections in the 4.3 deg 2 of A-ray coverage from Chan- 
dra and/or XMM telescopes (out of our full 6.7 deg 2 ). 
However, neither technique can exclude heavily-obscured 
AGNs effectively. 

Finally, IR colors can be used to exclude strong AGNs. 
In the IRAC color-color diagram with [3.6] — [4.5] versus 
[5.8] — [8.0], AGNs occupy a different locus from nor- 
mal galaxies, eit her SF or quiescent (jLacv et al.l I2004t 
IStern et al.ll2005[ ) . This is because the SED of an AGN in 
the IR can be well-represented by a power law, while that 
of a quiescent galaxy approximately follows Rayleigh- 
Jeans law and that of a SF galaxy is dominated by warm 
dust and PAH emission. To identi fy AGNs, w e adop t 
the empirical criteria suggested by IStern et all (|2005| ). 
These criteria may omit AGNs at high redshift, but are 
adequate for our purpose since our sample is at z < 0.5. 
In Figure [3j we show the IRAC color-color diagram for 
all the detections in both the 5.8 and 8.0 channels in 
the flu x-limited red-sequence sample, with IStern et al.l 
( 2005) criteria as the dotted wedge. Objects within the 
wedge are classified as AGNs (green filled squares). We 
also show SF and quiescent galaxies among the non- 
AGNs, defined in the next section, with blue filled tri- 
angles and red filled stars, respectively. The green dot- 
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Fig. 5. — Spectral templates shifted to z = 0.3, in AB magni- 
tudes, normalized to have the same magnitude zeropoint in the 
3.6 fim channel. W e use the template of IRAS12112 4- 0305 from 
IRieke et alj fl2009T ) to represent strong star-forming galaxies. For 
early- type galaxies without star formation and dust, we use a sim- 
ple stellar population (SSP) with an age of 10 Gyr and met allicity 
of O.4Z (SSP10G04Z) from the IBruzual fc Charloti 1120031 mod- 
els. This template is approximately the average of a variety of SSPs 
with age ranging from 5 to 10 Gyr and metallicity of OAZq and 
I.OZq. We show the average magnitudes in IRAC channels of these 
SSPs in black filled squares. To classify star-forming and quiescent 
galaxies, we define a demarcation template (Mixed) by linearly com- 
bining the two templates. Normalized at 2 (im, the Mixed template 
is the sum of 10% of IRAS12112 + 0305 and 90% of SSP10G04Z, 
and represents a galaxy with SSFR ~ 10 — 10 yr — 1 . The green dot- 
dashed line represents a power-law spectrum f ly) = u~ a of an 
AGN with a = 0.5. We also show the response curves of IRAC 
channels on the bottom. 

dashed line represents the colors of AGNs with power-law 
spectr a f(v) = v~ a ■, with a ranging from 0.5 to 3.0. 

The IStern et al.l ([2005D method to identify AGNs re- 
quires detections in all four channels. In the flux-limited 
red-sequence sample, a large fraction (1288 of 3310) are 
not detected in the 8.0 /im channel. Nevertheless, we will 
show in the next section that we can safely classify them 
as quiescent galaxies. Among those (2022 of 3310) de- 
tected in the 8.0 /an channel, 154 do not have matched 
detections in the 5.8 /im channel due to the relatively 
bright flux limits. For these objects, we use the flux limits 
of [5.8] < 19.8 (21.3) mag in the SWIRE (S-COSMOS) 
fields to calculate the lower limits of their [5.8] — [8.0] 
colors. We show the results in Figure |4] The colors of 
all 5.8 /im non-detections are unambiguously outside the 
AGN wedge and we therefore can safely classify them as 
non-AGNs and perform further SF and quiescent galaxy 
classification in the next section. 

3.4. Classification of star-forming and quiescent 
galaxies 

To separate SF and quiescent galaxies, we combine 
representative SED templates of strong SF galaxies and 
typical quiescent early- type galaxies. For strong SF 
galaxies, we choose the template of IRAS12112 + 0305 
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Fig. 6. — IR Color [3.6] — [8.0] of the spectral templates in Figure 
\E\ as a function of redshift. The various lines represent the colors 
of all the spectral templates in the same colors and line styles as 
in Figure [S] For comparison, we also show all the 8.0 fxra detec- 
tions in the blue-cloud (gray scale) and red-sequence galaxies (red 
dots) with i < 22.5 at z < 0.7. We classify galaxies above the 
Mixed template (the brown dashed line) as star-forming galaxies, 
if they are not in the AGN wedge in the IRAC color-color diagram 
(Figure El, At z > 0.5, the 6.2 lira PAH feature shifts out of the 
8.0 (im channel, causing the colors of star-forming and quiescent 
galaxies to blend into each other. 



from iRieke et al.l (|2009l ). an ultra-luminous IR galaxy 
(ULIRG) with sign ificant ongoing star formation (e.g., 
lArmus et a l. 2007). For typical quiescent early-type 
galaxies, we use a simple stellar population (SSP) with 
an age of 10 Gyr and m etallicity of OAZq from the 
iBruzual fc CharlotJ (|2003f ) models, with the Padoya, 1994 
library of stellar evolution tracks and the lChab ricr (2003) 
Initial Mass Function (IMF). We denote this template 
as SSP10G04Z. In Figure we show these templates, 
shifted to z = 0.3, within the wavelength range 2 to 
12 /mi. We also show a power-law spectrum f(v) = v~ a 
with a — 0.5 to represent a typical AGN. Note we 
have normalized all the SEDs to have the same mag- 
nitude zeropoint in the 3.6 fiiri channel for clarity. The 
IRAS12112 + 0305 template is dominated by warm dust 
continuum and strong PAH features. Meanwhile, the 
SSP template roughly follows the Rayleigh- Jeans law 
/(A) oc A~ 4 . This template is similar to the average 
of a variety of SSPs with age ranging from 5 Gyr to 15 
Gyr and metallicity of 0.4Z Q and 1.0Z Q . We show the 
average magnitudes of these SSPs in the IRAC channels 
in black filled squares. 

To define a demarcation template separating SF and 
quiescent galaxies, we linearly combine the templates of 
IRAS12112 + 0305 and SSP10G04Z. We first normalize 
them to have the same flux at 2 /xm (roughly if -band), 
and then define the demarcation mixed template (Mixed 
hereafter) to be the sum of 10% of IRAS12112 + 0305 
and 90% of SSP10G04Z. This Mixed template represents 
a galaxy with SSFR ~ 10~ 10 yr" 1 . We show the template 
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Fig. 7.— IR Color [3.6] - [8.0] vs. z for all 8.0 /im detections 
in the flux-limited sample (i < 20 mag in the SWIRE fields and 
i < 21 mag in the S-COSMOS field). Top panel: red-sequence 
galaxies. Many are scattered into the star-forming region (above 
the Mixed template). Middle panel: green-valley galaxies. Most of 
them appear in the star-forming region. Bottom panel: blue-cloud 
galaxies. Almost all of them are in the star-forming region. 

with the brown dashed line in Figure [S] 

Figure |5] presents the [3.6] — [8.0] colors of these tem- 
plates as a function of redshift and compares them with 
those of the 8.0 /im detections among the red-sequence 
(red dots) and blue-cloud galaxies (gray scale) with 
i < 22.5 at z < 0.7. The IRAS12112 + 0305 template 
clearly defines the upper edge of the space of the blue- 
cloud galaxies. Meanwhile, the Mixed template follows 
the lower edge of the blue-cloud galaxies. Most of the 
red-sequence galaxies have redder colors than the typical 
SSP template, SSP10G04Z, which should not be surpris- 
ing because very few galaxies are dust-free and can be 
accurately represented by SSPs. There are also many 
red-sequence galaxies scattered above the Mixed tem- 
plate, indicating they are in fact actively forming stars 
but their optical colors are reddened due to dust extinc- 
tion. We also see that above z = 0.5, the colors of SF 
and quiescent galaxies start to blend into each other be- 
cause the 6.2 /Ltm PAH emission band shifts out of the 
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Fig. 8. — IR Color [3.6] — [8.0] vs. z for non-detections in the 
8.0 (im channel in the flux-limited red-sequence sample. We use 
the flux limits [8.0] < 19.9 (21.0) mag in the SWIRE (S-COSMOS) 
fields to calculate the upper limits of their [3.6] — [8.0] colors. Al- 
most all the non-detections fall below the demarcation line (the 
Mixed template). There are four galaxies that can be identified as 
star-forming galaxies if their [8.0] fluxes are actually close to the 
flux limits. We treat them as quiescent galaxies, but identifying 
them as star-forming galaxies does not change our results. 

8.0 /im channel and we can not use the [3.6] — [8.0] color 
to identify SF galaxies any more. Also there are very few 
galaxies below z = 0.1 because of the small volume in the 
observation cone. We therefore only focus our analysis 
on the redshift range 0.1 < z < 0.5. 

In Figure [7J we present the [3.6] — [8.0] colors of all 
8.0 fim detections in the flux-limited sample split by op- 
tical color. In the top panel, we show that a significant 
fraction of the red-sequence galaxies are scattered into 
the SF region. In the middle panel, we show that most 
of the green- valley galaxies are actually actively forming 
stars at roughly the same level as the blue-cloud galaxies. 
Therefore, that the green valley is redder than the blue 
cloud is likely not an indication of a much lower SSFR, 
but may be due to dust reddening. In the bottom panel, 
we show almost all of the blue-cloud galaxies appear in 
the SF region. 

The classification of SF and quiescent galaxies requires 
detections in both the 3.6 /zm and 8.0 /nm channels. 
Among the 3310 red-sequence galaxies in the flux-limited 
sample, 1288 of them do not have matched 8.0 /im de- 
tections. For these galaxies, we use the flux limits of 
[8.0] < 19.9 (21.0) mag in the SWIRE (S-COSMOS) 
fields, to calculate the upper limits of their [3.6] — [8.0] 
colors. We show the results in Figure |8l Only four of 
these objects may be identified as SF galaxies if their 
[8.0] fluxes are close to the flux limits. Therefore the 
flux cuts in i band we adopt in Section 13.11 have pro- 
vided a sample with bright 3.6 /xm band fluxes (Figure 
[T]) and enable us to reliably classify almost all of the 
8.0 fim non-detections as quiescent galaxies. Below we 
treat the ambiguous four objects as quiescent galaxies 



Fig. 9.— Mo.: 



, h vs. z of the flux-limited red-sequence 



sample. The dotted line shows the flux cut i = 20 mag in the 
SWIRE fields and the dashed line represents the cut i = 21 mag 
in the S-COSMOS field, corrected with average iC-corrections of 
red-sequence galaxies. Due to Malmquist bias, we are investigating 
higher luminosity at higher redshift. 
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Fig. 10. — Fraction of obscured SF galaxies on the red sequence 
as a function of luminosity and redshift. In each redshift bin, we di- 
vide the sample into two subsamples with equal number of galaxies 
based on their luminosities and calculate the fraction of SF galax- 
ies and Poisson errors. At all redshifts, the fraction is higher at 
lower luminosity. We visualize the fraction with gray scale (darker 
color represents higher fraction). 
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TABLE 1 

Fraction of obscured star-forming galaxies on the red sequence 



z 


(Mo.3j - 51og 10 h) 


(Mo.3; - 51og 10 /l)cor a 


N KS b 


7V S F C 


Fraction of SF j 


0.1 < z < 0.2 


-18.47 


-18.64 


302 


99 


33 ± 4% 


(z) = 0.16 


-19.83 


-20.01 


304 


58 


19 ± 3% 


0.2 < z < 0.3 


-19.61 


-19.67 


504 


87 


17 ± 3% 


{2} = 0.25 


-20.77 


-20.83 


506 


50 


10 ± 2% 


0.3 < 2 < 0.4 


-20.20 


-20.15 


566 


91 


16 ± 2% 


{2} = 0.35 


-21.22 


-21.16 


567 


11 


8 ±2% 


0.4 < 2 < 0.5 


-20.83 


-20.65 


280 


50 


18 ± 3% 


(2} = 0.44 


-21.69 


-21.52 


281 


27 


10 ± 2% 



a (Mo.3 i — 51og 10 h) corrected to z — 0.3, assuming dMo.3 ■ / dz — —1.2 mag per unit redshift. 
k Number o£ red-sequence (RS) galaxies. 

c Number of star-forming (SF) galaxies in the red-sequence subsamplc. 



as well. However, we emphasize that treating them as 
SF galaxies does not have any noticeable effect in our 
analysis. We present the results in the next section. 

4. RESULTS: OBSCURED STAR FORMATION ON THE RED 
SEQUENCE 

With the classification scheme of AGNs, SF galaxies, 
and quiescent galaxies defined in Section [3j we present a 
quantitative study of the fraction of obscured SF galaxies 
on the red sequence in this section. 

4.1. Overall contribution 

In total, we have 3310 red-sequence galaxies in the 
flux-limited sample at 0.1 < z < 0.5. Among these we 
identify 15% (506) as SF galaxies, 84% (2783) as quies- 
cent galaxies, and 1% (21) as AGNs. Note if we did not 
exclude X-ray detections in about two-thirds of the to- 
tal 6.7 deg 2 fields, we would have 31 AGNs instead. In 
this flux-limited red-sequence sample, we find that over- 
all 15% have IR colors consistent with star formation, 
indicating their optical colors are likely reddened due to 
dust extinction. 

4.2. Luminosity dependence and redshift evolution 

It is well-known that elliptical galaxies dominate the 
bright end of the red sequence whil e disk-dominated sys- 
tems dominate the faint end (e.g., iMarinoni et al.|[l999l : 
iBundv et all201tt iZhu et alJl2010h . If the obscuration of 
the SF galaxies on the red sequence is due to high incli- 
nation, it is likely that the fraction of obscured SF galax- 
ies also depends on luminosity. Additionally, the global 
star formation rate (SFR) has decli ned by roughly an 
order-of-magnitude s ince z ~ 1 (e.g.. iMadau et alJll 996; 
Hopkins et al.l [20061 : IZhu et aT] I2009t IRu iopakar n et al.l 
2010T ). therefore we may expect that the fraction of ob- 



scured SF galaxies on the red sequence has declined with 
decreasing redshift. We investigate the luminosity de- 
pendence and redshift evolution in this section. 

For a flux-limited sample, due to Malmquist bias, the 
average luminosity is higher at higher redshift. We show 
this in Figure [SJ where we plot the absolute magnitude 
— 5 log 10 h versus z of the red-sequence galaxies. 
We also show the i band flux cuts, 20 mag in the SWIRE 
fields and 21 mag in the S-COSMOS field, corrected with 
average i^-corrections of red-sequence galaxies, with the 
dotted and dashed lines, respectively. 
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Fig. 11. — Fraction of obscured SF galaxies on the red sequence 
as a function of luminosity and redshift, as in Figure ITOl We show 
the fraction at the mean magnitude in each subsample. In the 
bottom panel, we assume that 0,3 i band absolute magnitude fades 
by 1.2 per unit redshift and correct all magnitudes to z = 0.3. In all 
redshift bins, the fraction (in percent) is higher at lower luminosity 
by ~ 8% per mag. Assuming linear luminosity dependence, at 
L > L* , the percentage of obscured SF galaxies has decreased by 
~ 6 ± 3% from z = 0.44 to z = 0.25. At L < L* , the data from 
2 = 0.35 to 2 = 0.16 are consistent with each other, within the 
errors. At fainter luminosity (> 1 mag fainter than L*), we do not 
have a wide enough redshift baseline to test for redshift evolution. 
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To overcome the Malmquist bias and possible degen- 
eracy between luminosity and redshift dependence, we 
therefore define four redshift bins with a small binsize 
Az = 0.1. In each redshift bin, we divide the sample 
into two subsamples with equal number of red-sequence 
galaxies based on their — 51og 10 h magnitudes. 

Figure [TU] presents the fraction of SF galaxies with 
Poisson errors at the mean absolute magnitude and red- 
shift of each subsample. We visualize the fraction with 
gray scale (darker colors indicate higher fractions). We 
also list the results in Table [U The fraction appears to 
be higher at lower luminosity in all four redshift bins. 

We show the luminosity dependence in Figure 111! 
where we plot the fraction as a function of the mean 
Mo.3i — 5 log 10 h in each subsample in the top panel. In 
the bottom panel, we assume that the luminosity fades 
with decreasing redshift by 1.2 mag per unit redshift, 
which is the lu minosity evolution factor for an SSP with 
Salpeterl (Il955h IMF formed at z ~ 4 (e.g., Ivan Dokkuml 
2008), and correct all luminosities to z — 0.3. We 
find that on average, at 0.1 < z < 0.5, the fraction 
of SF galaxies on the red sequence is about 15% at 
Mo.si - 51og 10 /i ~ -20.5 (about L13)- And at all red- 
shifts, the fraction (in percent) increases by ~ 8% per 
mag with decreasing luminosity to 0.2L*. For example, 
in the redshift range 0.3 < z < 0.4, the fraction increases 
from 8±2% to 16±2% as luminosity decreases from 
Mo.3i - 51og 10 /i 21.2 to -20.2. 

Also due to Malmquist bias, at different lumi- 
nosities we a re iny estigating different redshift ranges. 
iBundv et all (|2010l) find that at stellar mass < lO n M , 
the abundance of red-sequence disk galaxies increases 
with decreasing redshift, but not as fast as red-sequence 
spheroidals. They also find that at higher stellar mass, 
the disk population has declined since z ~ 1. If the ob- 
scured SF population behaves as their selected disk pop- 
ulation, we expect that the redshift evolution of the frac- 
tion also depends on luminosity. With our flux-limited 
sample, each luminosity range covers a different redshift 
range, complicating our attempt to investigate this ques- 
tion. Figure[TT]does not reveal any strong trend of evolu- 
tion strength with luminosity, but we cannot address this 
question substantially below L* with this sample. With 
our sample, assuming linear luminosity dependence, at 
L L* (e.g., at — b\og w h ~ —21), the percent- 

age of obscured SF galaxies on the red sequence has de- 
creased by ~ 6 ± 3% from z = 0.44 to z = 0.25. At 
L < L*, the data from z = 0.35 to z = 0.16 are consis- 
tent with each other, within the errors. At fainter lumi- 
nosity (> 1 mag fainter than L*), we do not have a wide 
enough redshift baseline to put interesting constraints on 
evolution. 

4.3. The origin of 24 fim emission on the red sequence 

We have shown that a significant fraction (15%) of red- 
sequence galaxies have IR colors consistent with star for- 
mation. Another approach to search for star formation 
is to use the 24 /jm band of MIPS, also onboard Spitzer. 
The 24 /x m emission is a close tracer of SFR of S F galax- 
ies (e.g., iCalzetti et alj|2007t iRieke et all 120091). Previ- 
ous studies (e.g.. iDavoodi et al.ll2006HRo dighicro et al.l 

13 At z ~ 0.1, the L* is M .i r - 51og 10 /i ~ -20.4 in 01 r band 
(Blanton ct al. 2003b) and 01 r band is very close to 3 i band. 
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Fig. 12. — IR color [3.6] — [8.0] vs. z for red-sequence galaxies with 
24 /im emission f24 > 300 fijy. We show star-forming galaxies in 
blue filled triangles, quiescent galaxies in red filled stars, and AGNs 
in green filled squares. At 0.1 < z < 0.5, there are 355 (~ 11% of 
the sample) galaxies with f'24 > 300 fiJy, among which wc identify 
~ 74% (265) as SF galaxies, ~ 23% (81) as quiescent galaxies, and 
~ 3% (9) as AGNs. The various lines represent spectral templates, 
as in Figure [5] 



[2CMlBrand"et al. 2009) find - 10-20% of red-sequence 
galaxies with 24 /j,m emission /24 ^ 300 fiJy. With 
the 24 /zm imaging data, we also find that in our flux- 
limited red-sequence sample at 0.1 < z < 0.5, about 
11% (355) have 24 fj,m flux / 24 > 300 /zJy. At z = 0.3, 
this 24 ^m flux corresponds to a SFR ~ 3 M yr _1 
(|Rieke et al.|[2009h . 

The 24 /x m emission, however, can also originate from 
AGNs (e.g., iFarrah et aLll2003t Brand et al.l l2006). It is 
interesting to investigate whether star formation or AGN 
or both activities are responsible for the 24 /im emission 
of these red-sequence galaxies. Previous studies (e.g., 
Davoodi et "aT1l2006HRodighiero et aT]|2007l: lB7and et al.l 



20091 ) find that both star f ormation and AGN activ- 
ities could be responsible. iBrand et al.l (|2009l ) select 
89 (~ 10%) red-sequence galaxies with / 2 4 > 300 /iJy 
from the AGN and Galaxy Evolution Survey (AGES) 
and employ both th e emission-line diagnostic diagram 
(|Baldwin et al.lll981l BPT) and the IRAC color-color di- 
agram to investigate the origin of the emission. They find 
that in their sample, with the former method the major- 
ity (~ 64%) are identified as AGNs, while with the latter 
they identify ~ 64% as SF galaxies, ~ 31% as quiescent 
galaxies, and only ~ 5% as AGNs. Therefore both star 
formation and AGN activities take place in their sam- 
ple, but in the optical the star formation is obscured by 
dust and in the IR the AGN emission is outshone by star 
formation emission. 

The PRIMUS spectroscopy does not resolve the narrow 
lines required to construct the BPT diagram. However, 
we can perform the same AGN, SF galaxy, and quiescent 
galaxy classification as above using IRAC colors. Among 



10 



Zhu et al. 



the 355 red-sequence galaxies with /24 > 300 //Jy, we 
identify ~ 74% (265) as SF galaxies, ~ 23% (81) as qui- 
escent galaxies, and ~ 3% (9) as AGNs. In Figure [T2l we 
plot the [3.6] — [8.0] colors of these galaxies as a function 
of red shift. Our results are consistent with IBrand et al] 
(2009) and suggest that star formation should account 
for at least part of the 24 /j,m emission. 

5. DISCUSSION 

5.1. The nature of galaxies on the red sequence 

Although the red sequence is a well-defined locus with 
very small scatter in the color-magnitude diagram, it ac- 
tually consists of both early-type and late-type galaxies: 
early-type spheroidal galaxies, late-type disk-dominated 
galaxies without ongoing star formation, and late-type 
dusty SF galaxies with their optical colors reddened due 
to dust extinction (e.g., edge-on Sb/Sc, etc.). When 
studying galaxy evolution based on optical color, it is 
therefore necessary to understand the contribution of dif- 
ferent types of galaxies on the red sequence. Here we 
describe a number of methods of doing so, finding that 
broadly speaking they agree with our results here on the 
fraction of star- forming and/or disk galaxies on the red 
sequence. 

There are many ways to study the composition of 
galaxies on the red sequence. A non-exhaustive list in- 
cludes: 



(1) 



direct c lassification into diff e rent morphologies 
(e.g.. IMarinoni et alJ 119991: lYan fc Thompson 

120041 IBell et al 
lLotz et all 1 20081: 
iBundv et ail bold 




Moustakas et al. 
iWeiner et all l200l 
Blanton fc Moustakas! 2009; 



Zhu et aLH20Trf ~ 



(2) 



using inclination-corrected properties (e.g., 
Tully et all 119981: [Masters et all 120031: iShao et all 
20071: iDriver et all 120071: IBailin fc Harris! l2008aBrI 
Unterborn fc Rvden 2008; Maile r et alll2009D ; 



(3) SED fitti ng where dust extinction is a free parame- 
ter (e.g.. IKauffmann et al J 120031: 1 Wolf et al.l l200l 
iBrammer et al.ll2009l ): 

(4) using optical color and o ptical-near-IR (NIR) 
color-c olor diagram (e.g., iPozzetti fc Mannuccil 
20001: lLabbe et all 12003 : iWuvts et all 120071 : 
Williams et al.ll2009l ): 

IBrand et al1l2009l this work) 

JO 



(5) using IR color (e.g. 



(6) 



(7) 



using MIR 



Davoodi et al 
Saintonge et al.l 




2005 



work) ; 
using 



e.g^ moia et all 

iRodighiero et al| _ 
IBrand et all 120091 this 



2007 



emission line measurements 



Hammer et al. 


1997; Dcmarco et all 


Popesso et all 


20071: IVerdugo et al. 20081 



2005; 



It is difficult to carry out an exact comparison of these 
studies because of different sample selection methods, 
different luminosity and mass ranges investigated, and 
different redshift ranges investigated. However, all of 
these studies lead to the conclusion that a significant 
fraction of red-sequence galaxies are disk galaxies and/or 
obscured SF galaxies. Detailed studies conclude that the 
fraction depends strongly on luminosity and/or mass. In 



terms of morphology, in the local universe, va rious papers 
(e.g., IMarinoni et allll999l : IZhu et alll2010l ) have shown 
that spheroidal galaxies dominate the bright end (> L*) 
and disk-dominated systems d ominate at the faint end. 
At redshifts 0.2 < z < 1.2, IBundv et all (poTTl find 
that disk galaxies represent nearly one-half of all red- 
sequence population and dominate at low stellar mass 
< 10 10 Mq. In terms of star formation, we find at red- 
shifts 0.1 < z < 0.5, the fraction of obscured SF galaxies 
on the red sequence increases by ~ 8% per mag with 
decreasing luminosity to L ~ 0.2L*. 

The redshift evolution of the relative fractions of dif- 
ferent comp_oneirts_ojithfi_red sequen ce, however, is less 
clear Mousta kas et all (|2004j . see also lYan fc Thompson! 
2003) find that the fraction of spheroidal galaxies in ex- 
tremely red objects (EROs) is 33% - 44% at redshift 
~ 1-2. lLotz et all (|2008fl find that at luminosity 
brighter than 0AL* B , edge-on and dusty disk galaxies 
are almost one third of the red sequence at z ~ 1.1, 
while they only make up < 10% at z ~ 0.3. We also 
find that at L > L*, the fraction of obscured SF galax- 
ies decreases by ~ 6 ± 3% from z — 0.44 to z = 0.25. 
However, ot her studies do not find a noticeable redshift 
dependence. IBell et all (|2004b[ ) and IWeiner et al.l (|2005f l 
find that the fraction of edge-on disk galaxies is similar at 
z ~ 1 and at z ~ 0. One possible reason for the disagree- 
ment is that the redshift evolution of relative fractions 
may depend on other pr operties such as ste llar mass or 
luminosity. For example, IBundv et all ([20101 ) find that at 
stellar mass < 1O 11 M , the abundance of red-sequence 
disk galaxies increases with decreasing redshift, but not 
as fast as red-sequence spheroidals; meanwhile, at higher 
stellar mass, the disk population has declined since z ~ 1. 

The relative fractions of different types of galaxies and 
their redshift evolution on the red sequence may also de- 
pend on environment. It is well-known that overall early- 
type galaxies are more frequently found in high-density 
environments while late-type galaxies tend to dominate 
in low-density regions (the Morphology-density relation, 
lDresslerlll9 80). There is also the SSFR-density relation 
such that the mean SSFR is higher at low density, an ef- 
fect that is impo rtant in the local universe as well as 
at z ~ 1 (e.g., IKauffmann et al.l 120041 : iCooper et all 
2008). Meanwhile, the mean SFR is higher at low den- 
sity in the local universe but the SFR-den sity relation 
appe a rs to be reversed a t z ~ 1 (e.g., lElbaz et"al"1 
120071: ICooper et all I2008D . The fraction of obscured 
SF gal axies therefore may also depend on environment 
(e.g.. IWolf et all [20051: IGallazzi et all [20081 : IWolf et all 
120091: iFeruglio et all 120101 ) . These studies suggest that 
in interpreting the optical color-magnitude diagram as a 
function of environment, the population of obscured SF 
galaxies may need to be accounted for. 

An interesting suggestion by some papers is that 
the obscured SF galaxies on the red sequence are 
a third type ot her than red and blue galaxies (e.g.. 
Wolf et all 120051: iWeinmann et al l 120061: LPopesso et all 
20071: IBailin fc Harris! I2008at IVerdugo et all [20081: 
Wolf et al l I 2009D. For example. IWolf et all (|2005l ) and 
Wolf et all (120091 ) find these galaxies are distributed 



differently in environment than both the red and the 
blue galaxy populations: unlike red galaxies they tend 
to avoid dense cluster regions, but in contrast to blue 
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Fig. 13. — Color-magnitude diagram for galaxies between z = 0.2 and z ■■ 
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but only for red-sequence galaxies. The bimodality is more striking using IR color. And a significant fraction of red-sequence galaxies are 
scattered into the star-forming region in the IR color-magnitude diagram. 



galaxies they are not common in low- density fields. 
With 24 fim MIR data. IWolf eTail (p009t ) further show 
that the mean SSFR of these galaxies in clusters is lower 
than in the field, in contrast to that of blue galaxies 
alone, which appears similar in different environments. 
They also show using HST imaging that edge-on spirals 
form only a small fraction of spiral galaxies. These 
authors thus suggest that these obscured SF galaxies are 
a transitioning population between blue field galaxies 
and red cluster galaxies in cores. 

Finally, an interesting question is how much these ob- 
scured SF galaxies a ccount for in the overall SFR density 
in the universe (e.g., iBell et al.ll2007l ). In the flux-limited 
sample, we have 506 obscured SF galaxies on the red se- 
quence, 1451 green-valley and 2497 blue-cloud galaxies. 
Assuming that they have comparable average SSFR and 
stellar mass, we roughly estimate that of order 10% of to- 
tal cosmic star formation is contributed by obscured SF 
galaxies. An exact analysis requires precise SFR mea- 
surements, beyond the scope of this paper. We will ad- 
dress this issue in future papers measuring the global 
SFR density using UV and 24 /im MIR luminosity. 

5.2. IR color-magnitude diagram 

As a proxy for the SSFR-stellar mass relation, the opti- 
cal color-magnitude diagram is a powerful tool in study- 
ing galaxy evolution; however, we argue here that IR 
color can be more effective in separating star-forming and 
quiescent galaxies when it is available. As shown above, 
optical color can be a misleading variable for SSFR and 
there is a significant fraction of obscured SF galaxies on 



the red sequence. The fraction depends on luminosity 
and redshift and may also depend on environment. This 
complication affects galaxy evolution studies based on 
optical color. Therefore, it is useful to find more suit- 
able variables tracing SSFR. Exampl es of such variables 
are inclination-corrected color (e.g., M ailer et al.ll2009T) 
and extinction-correc ted color from SED fitting (e.g., 
iBrammer et al.l [2009( 1 . IR color is another, and better, 
variable than optical color to represent SSFR; therefore, 
replacing optical color with IR color when available is a 
better strategy when constructing the color-magnitude 
diagram. 

As an example, in Figure 13 we show both the opti- 
cal color-magnitude diagram (upper-left panel) and IR 
color- magnitude diagram (lower-left panel), where we 
have reversed the IR color [3.6] - [8.0] to [8.0] - [3.6] 
so that SF galaxies are on the bottom. We show all the 
galaxies in the flux-limited sample at 0.2 < z < 0.4. 
For non-detections in 8.0 /^m channel, we use the lower 
limits by applying the flux limits [8.0] < 19.9 (21.0) 
mag in the SWIRE (S-COSMOS) fields. We have also 
applied simple i^-corrections using the templates for 
IRAS12112 + 0305 and SSP10G04Z to correct the IR 
color [8.0] — [3.6] so that it is the color in the channels 



shifted blueward to 



0.3, 



Mo, 



[8,0] 



Mo, 



[3.6] - 



The bimodality in the IR color-magnitude diagram is 
more striking than in the optical color-magnitude di- 
agram. This should not be surprising because the bi- 
modality is caused by the short time-scale of transition 
from star formation to quiescence. In the right panels, 
we only plot the red-sequence galaxies and show that a 
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significant fraction of them are in the SF region in the 
IR color-magnitude diagram, as we quantified in detail 
above. 



6. SUMMARY 

With the recently completed PRIsm MUlti-object Sur- 
vey (PRIMUS), we perform a quantitative study of the 
fraction of obscured star-forming galaxies on the red se- 
quence at intermediate redshift. PRIMUS targeted fields 
with existing deep infrared (IR) imaging from the IRAC 
and MIPS instruments onboard the Spitzer Space Tele- 
scope, including 5.7 deg 2 fields from the SWIRE survey 
and a 1 deg 2 field from the S-COSMOS survey. With 
the precise redshifts from PRIMUS, we select in these 
fields an i band flux-limited sample of 3310 red-sequence 
galaxies, with i < 20 mag in the SWIRE fields and i < 21 
mag in the S-COSMOS field, at 0.1 < z < 0.5. With the 
IR imaging data, we classify the red-sequence galaxies 
into AGNs, star-forming galaxies, and quiescent galax- 
ies. Our sample is the largest sample at intermediate 
redshift for such a study and we present for the first 
time a quantitative analysis of the fraction of obscured 
star-forming galaxies on the red sequence as a function 
of luminosity. 

We find that on average, at L ~ L* , about 15% of 
red-sequence galaxies have IR colors ([3.6] — [8.0]) con- 
sistent with star-forming galaxies at 0.1 < z < 0.5. The 
percentage of obscured star- forming galaxies increases by 
~ 8% per mag with decreasing luminosity to L ~ 0.2L*. 
At L > L* , the fraction has declined by ~ 6 ± 3% from 
z = 0.44 to z = 0.25. At L < L* , the fraction is con- 
sistent with no redshift evolution between z = 0.35 and 
z = 0.16. 

Using the 24 ^m imaging from MIPS on board Spitzer, 
we find that ~ 11% of red-sequence galaxies exhibit 
24 ijluy emission with / 2 4 > 300 /xJy at redshifts 0.1 < 
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z < 0.5, in agreement with previous studies. With the 
IRAC imaging, we find that ~ 74% of these galaxies have 
IR colors consistent with star-forming galaxies, suggest- 
ing that star formation accounts for at least part of the 
24 fim. emission. 

Our analysis suggests that a significant fraction of red- 
sequence galaxies are actually actively forming stars, but 
their optical colors are reddened due to dust extinction. 
An accurate understanding of galaxy evolution using op- 
tical color therefore needs to account for this complica- 
tion. We suggest that IR color, e.g., the IRAC color 
[3.6] — [8.0], should be a better color to use when con- 
structing the color-magnitude diagram to represent the 
specific star-formation rate-stellar mass relation. 

We acknowledge Rebecca Bernstein, Adam Bolton, 
Douglas Finkbeiner, David W. Hogg, Timothy McKay, 
Sam Roweis, and Wiphu Rujopakarn for their contri- 
butions to the PRIMUS project. This paper includes 
data gathered with the 6.5 meter Magellan Telescopes 
located at Las Campanas Observatory, Chile. We thank 
the support staff at LCO for their help during our 
observations, and we acknowledge the use of commu- 
nity access through NOAO observing time. Funding 
for PRIMUS has been provided by NSF grants AST- 
0607701, 0908246, 0908442, 0908354 and NASA grant 
08-ADP08-0019. RJC is supported by NASA through 
Hubble Fellowship grant HF-01217 awarded by the Space 
Telescope Science Institute, which is operated by the As- 
sociated of Universities for Research in Astronomy, Inc., 
for NASA, under contract NAS 5-26555. 

This work is based in part on observations made with 
the Spitzer Space Telescope, which is operated by the 
Jet Propulsion Laboratory (JPL), California Institute of 
Technology under a contract with NASA. We would like 
to thank the SWIRE and the S-COSMOS teams for mak- 
ing their data publicly available. 



Cimatti, A., Daddi, E., & Renzini, A. 2006, A&A, 453, L29 
Coia, D., et al. 2005, A&A, 430, 59 
Colless, M., et al. 2001, MNRAS, 328, 1039 
Coil, A. L., et al. 2010, ApJ, submitted 
Cool, R. J., et al. 2008, ApJ, 682, 919 
Cooper, M. C, et al. 2008, MNRAS, 383, 1058 
Davis, M., Huchra, J., Latham, D. W., & Tonry, J. 1982, ApJ, 
253, 423 

Davis, M., et al. 2003, in Presented at the Society of 

Photo-Optical Instrumentation Engineers (SPIE) Conference, 
Vol. 4834, Society of Photo-Optical Instrumentation Engineers 
(SPIE) Conference Series, ed. P. Guhathakurta, 161-172 

Davoodi, P., et al. 2006, MNRAS, 371, 1113 

Demarco, R., ct al. 2005, A&A, 432, 381 

Dressier, A. 1980, ApJ, 236, 351 

Driver, S. P., Popescu, C. C, Tuffs, R. J., Liske, J., Graham, 

A. W., Allen, P. D., & de Propris, R. 2007, MNRAS, 379, 1022 
Eisenhardt, P. R., et al. 2004, ApJS, 154, 48 
Elbaz, D., et al. 2007, A&A, 468, 33 
Faber, S. M., et al. 2007, ApJ, 665, 265 

Farrah, D., Afonso, J., Efstathiou, A., Rowan-Robinson, M., Fox, 

M., & Clements, D. 2003, MNRAS, 343, 585 
Fazio, G. G., et al. 2004, ApJS, 154, 10 
Feruglio, C, et al. 2010, ApJ, 721, 607 

Gallazzi, A., Brinchmann, J., Chariot, S., & White, S. D. M. 

2008, MNRAS, 383, 1439 
Garilli, B., et al. 2008, A&A, 486, 683 
Hammer, F., et al. 1997, ApJ, 481, 49 



PRIMUS: Obscured Star Formation 



13 



Hopkins, P. F., Hcrnquist, L., Cox, T. J., Di Matteo, T., 
Robertson, B., & Springel, V. 2006, ApJS, 163, 1 

Kauffmann, G., White, S. D. M., Heckman, T. M., Menard, B., 
Brinchmann, J., Chariot, S., Tremonti, C, & Brinkmann, J. 
2004, MNRAS, 353, 713 

Kauffmann, C, et al. 2003, MNRAS, 341, 33 

Labbe, I., et al. 2005, ApJ, 624, L81 

Lacy, M., et al. 2004, ApJS, 154, 166 

Li, A., & Draine, B. T. 2001, ApJ, 554, 778 

Lilly, S. J., et al. 2007, ApJS, 172, 70 

Lonsdale, C. J., et al. 2003, PASP, 115, 897 

Lotz, J. M., et al. 2008, ApJ, 672, 177 

Madau, P., Ferguson, H. C, Dickinson, M. E., Giavalisco, M., 
Steidel, C. C, & Fruchter, A. 1996, MNRAS, 283, 1388 

Mailer, A. H., Berlind, A. A., Blanton, M. R., & Hogg, D. W. 
2009, ApJ, 691, 394 

Marinoni, C, Monaco, P., Giuricin, G., & Costantini, B. 1999, 
ApJ, 521, 50 

Martin, D. C, ct al. 2005, ApJ, 619, LI 

Masters, K. L., Giovanelli, R., & Haynes, M. P. 2003, AJ, 126, 158 

Moustakas, L. A., et al. 2004, ApJ, 600, L131 

Popesso, P., Biviano, A., Romanicllo, M., & Bohringcr, H. 2007, 

A&A, 461, 411 
Pozzetti, L., & Mannucci, F. 2000, MNRAS, 317, L17 
Rieke, G. H., Alonso-Herrero, A., Weiner, B. J., Perez-Gonzalez, 

P. G., Blaylock, M., Donley, J. L., & Marcillac, D. 2009, ApJ, 

692, 556 

Rieke, G. H., et al. 2004, ApJS, 154, 25 
Rix, H., et al. 2004, ApJS, 152, 163 

Rodighicro, G., Cimatti, A., Franceschini, A., Brusa, M., Fritz, 
J., & Bolzonella, M. 2007, A&A, 470, 21 



Rujopakarn, W., et al. 2010, ApJ, 718, 1171 

Saintonge, A., Tran, K., & Holden, B. P. 2008, ApJ, 685, L113 

Salpeter, E. E. 1955, ApJ, 121, 161 

Sanders, D. B., et al. 2007, ApJS, 172, 86 

Shao, Z., Xiao, Q., Shen, S., Mo, H. J., Xia, X., & Deng, Z. 2007, 

ApJ, 659, 1159 
Smith, J. D. T., et al. 2007, ApJ, 656, 770 
Stern, D., et al. 2005, ApJ, 631, 163 
Strateva, I., et al. 2001, AJ, 122, 1861 

Tully, R. B., Pierce, M. J., Huang, J., Saunders, W., Verheijen, 

M. A. W., & Witchalls, P. L. 1998, AJ, 115, 2264 
Unterborn, C. T., & Ryden, B. S. 2008, ApJ, 687, 976 
van Dokkum, P. G. 2008, ApJ, 674, 29 

Verdugo, M., Ziegler, B. L., & Gerken, B. 2008, A&A, 486, 9 
Weiner, B. J., et al. 2005, ApJ, 620, 595 

Weinmann, S. M., van den Bosch, F. C, Yang, X., & Mo, H. J. 

2006, MNRAS, 366, 2 
Werner, M. W., et al. 2004, ApJS, 154, 1 

Williams, R. J., Quadri, R. F., Franx, M., van Dokkum, P., & 

Labbe, I. 2009, ApJ, 691, 1879 
Willmer, C. N. A., ct al. 2006, ApJ, 647, 853 

Wolf, C, Gray, M. E., & Meisenheimer, K. 2005, A&A, 443, 435 
Wolf, C, Meisenheimer, K., Rix, H., Borch, A., Dye, S., & 

Kleinheinrich, M. 2003, A&A, 401, 73 
Wolf, C, et al. 2004, A&A, 421, 913 
— . 2009, MNRAS, 393, 1302 
Wuyts, S., et al. 2007, ApJ, 655, 51 
Yan, L., & Thompson, D. 2003, ApJ, 586, 765 
York, D. G., et al. 2000, AJ, 120, 1579 

Zhu, G., Blanton, M. R., & Moustakas, J. 2010, ApJ, 722, 491 
Zhu, G., Moustakas, J., & Blanton, M. R. 2009, ApJ, 701, 86 



